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Establishing a connection between high-power pulsars and very-high-energy 
gamma-ray sources 
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Abstract: In the very-high-energy (VHE) gamma-ray wave band, pulsar wind nebulae (PWNe) represent 
to date the most populous class of Galactic sources. Nevertheless, the details of the energy conversion 
mechanisms in the vicinity of pulsars are not well understood, nor is it known which pulsars are able to 
drive PWNe and emit high-energy radiation. In this paper we present a systematic study of a connection 
between pulsars and VHE 7-ray sources based on a deep survey of the inner Galactic plane conducted 
with the High Energy Stereoscopic System (H.E.S.S.). We find clear evidence that pulsars with large 
spin-down energy flux are associated with VHE 7-ray sources. This implies that these pulsars emit on the 
order of 1% of their spin-down energy as TeV 7-rays. 



In 1989, the Crab Nebula was discovered as the 
first celestial source ofVHE 7-radiation [15]. The 
pulsar inside the nebula drives a powerful wind of 
highly relativistic particles that ends in a termina- 
tion shock from which high-energy particles with 
a wide spectrum of energies emerge [8]. High- 
energy electrons' among these particles can give 
rise to two components of electromagnetic radia- 
tion: a low-energy component from synchrotron 
radiation and a high-energy component from in- 
verse Compton (IC) up-scattering of ambient pho- 
tons. 

Recently, advances in VHE instrumentation have 
made the discovery of many new, predominantly 
Galactic, sources possible. Of these, a signifi- 
cant number can be identified as PWNe. Promi- 
nent examples are the PWN of the energetic pul- 
sar PSR B 1509— 58 in the supernova remnant 
MSH 15-52 [2], and HESS J0835-455 [3], as- 
sociated with Vela X, the nebula of the Vela pulsar. 

These 7-ray PWNe are extended objects with an 
angular size of a fraction of a degree, translating 
into a size of some lOpc for typical distances of a 
few kpc. In addition to the open puzzle of pulsar 



spin-down power conversion, a surprising obser- 
vation is that the centroids of these 7-ray PWNe 
are often displaced from their pulsars by distances 
similar to the nebular size. Such displacements, 
although usually at smaller scales, are also seen in 
some X-ray PWNe. The origin of the displacement 
remains unknown. It might be attributed to pulsar 
motion (e.g. [14]), causing the pulsar to leave its 
nebula behind, or to a density gradient in the am- 
bient medium [6]. 

The aforementioned examples of coincidences be- 
tween VHE 7-ray sources and radio pulsars mo- 
tivated a systematic search for VHE counterparts 
of energetic pulsars using the H.E.S.S. system of 
imaging Cherenkov telescopes located in Namibia 
[9]. To be detectable by H.E.S.S., a source 
at distance d has to provide a 7-ray luminos- 
ity in the 1 TeV to 10 TeV range of L-^ ~ 
lO'^^ (P erg s^'^kpc^^. Assuming a conversion ef- 
ficiency of 1 % of pulsar spin-down energy loss E 
into TeV 7-rays (where E is determined from the 
measurement of the rotation period ft and the rate 



1. here and in the following, 'electrons' refers to both 
electrons and positrons 
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at which the rotation slows down O), PWNe of pul- 
sars with E around 10^* (P ergs^^kpc"^ might 
be detectable. We note that for typical electron 
spectra, only a small fraction of the total energy 
in electrons is carried by the multi-TeV electrons, 
that are responsible for TeV 7-rays by IC scattering 
off ambient photons (including those from the cos- 
mic microwave background) and for ke V ^ -ray s by 
synchrotron radiation. Even a 1% energy output in 
TeV 7-rays already imphes a large fraction of spin- 
down energy loss going into relativistic electrons. 

Here we investigate how the probability to detect 

in VHE 7-rays PWNe surrounding known pulsars 
varies with the spin-down energy loss of the pul- 
sar, testing the plausible assumption that the 7-ray 
output of a PWN correlates in some fashion with 
the power of the pulsar feeding it. 

The VHE 7-ray data set used to search for 7-ray 
emission near the location of known radio pul- 
sars comprises all data used in the H.E.S.S. Galac- 
tic plane survey [1, 5], including an extension of 
the survey to Galactic longitudes —60° < I < 
—30°, dedicated observations of Galactic targets 
and re-observations of H.E.S.S. survey sources. 
The search covers a range in Galactic longitude 
from —60° to 30° while the range in Galactic lat- 
itude is restricted to ±2deg, a region well cov- 
ered in the survey. A total of 435 pulsar locations 
are tested, taken from the Parkes Multibeam Pul- 
sar Survey (PMPS, [10] and references therein), 
as recorded in the ATNF pulsar catalogue. Pul- 
sars without measured period derivatives are ig- 
nored. Over the range of the H.E.S.S. survey, 
the PMPS provides reasonably uniform sensitivity 
[11], enabhng a rehable estimate of the frequency 
of chance coincidences between a 7-ray source and 
a pulsar. The analysis of the 7-ray data follows the 
standard H.E.S.S. analysis [4]. Initially, a sky map 
is generated providing the significance of a 7-ray 
excess for a given position. Taking into account 
the properties of known 7-ray PWNe, the search 
is optimised for sUghtly extended sources - on the 
scale of the angular resolution (« 0. 1 deg) of the 
H.E.S.S. telescopes - and allows for small offsets 
from the pulsar positions. Each excess is deter- 
mined by counting 7-ray candidate events within 
e < 0.22 deg (6*2 < 0.05 deg^) of a given position 
and subtracting a background estimated from areas 
in the same field of view. The sky map is used to 



look up the significance of a 7-ray excess at the po- 
sition of the radio pulsars, as well as for randomly 
generated test positions used to evaluate the sta- 
tistical significance of the association (details are 
given below). We require an excess significance of 
at least 5 standard deviations above the background 
as a signature of a VHE 7-ray signal. Given the 
modest number of trials - the 435 pulsar locations 
- the number of false detections is negligible with 
this requirement and in any case small compared 
to the probability for chance coincidences between 
radio pulsars and VHE 7-ray sources. 
Of the 435 pulsars, 30 are found with signifi- 
cant 7-ray emission at the pulsar location (Fig. 1, 
top left panel). The lower left panel of Fig. 1 
displays the fraction of pulsars with such 7-ray 
emission for different intervals in spin-down flux 
E/(f. The fraction is about 5% for pulsars with 
spin-down flux below 10^^ ergs~^kpc~^ and in- 
creases to about 70% for pulsars with E/cP above 
10^^ ergs~^kpc~^. Not all of these associations 
are necessarily genuine. The rate of chance coinci- 
dences is estimated by generating 10^ realisations 
of random pulsar samples (each consisting on aver- 
age of 435 "pulsars") following the distribution in 
longitude and latitude of the PMPS pulsars and tak- 
ing into account the narrowing of the distribution 
in latitude with increasing spin-down flux. The ex- 
pected fraction of chance coincidences is shown 
as dark shaded areas in Fig. 1 and varies between 
4% to 12%. All associations with pulsars with 
E/(P < 10'^^ ergs^^kpc^^ are within statistical 
errors consistent with chance coincidences. Indeed 
for plausible values of the ratio between the 7-ray 
luminosity and the pulsar spin-down energy loss, 
L^/E, no detectable emission would be expected 
from such pulsars. On the other hand, the detection 
of emission from high-power pulsars is statistically 
significant. The probability that the detection of 
VHE sources coincident with 9 or more of the total 
of 23 pulsars above £^/c?^ > 10'^'' erg s^^kpc^^ re- 
sults from a statistical fluctuation is ^ 3.4 x 10~^. 
For detection of 5 or more of the total of 7 pulsars 
above 10^^ erg s~^kpc~^, the chance probabihty is 
~ 4.2 X 10"*. 

Given the high density of pulsars, a single 7-ray 
source may even coincide with more than a single 
pulsar, and thus appear more than once amongst 
the "detections" in the upper left panel of Fig. 1. 
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Figure 1: Top row: Distribution in logio(i?/(i'^) of all PMPS pulsars in the H.E.S.S. scan range (shaded in 
light grey), of chance coincidences (shaded in dark grey) and of detected pulsars (black line). Here, E/(P 
is measured in erg s^^kpc^^. Bottom row: The points show the fraction of pulsars with significant 7-ray 
excess at the pulsar position, as a function of logio(£'/(i^). The shaded band represents the probability 
for a chance coincidence. The width of the band accounts for the uncertainty in the width of the latitude 
distribution of pulsars. Left: all pulsars; right: double occurrences of gamma-ray sources removed by 
omitting pulsars which overlap with stronger pulsars or known non-pulsar sources. 



Removal of such double occurrences (Right panels 
of Fig. 1) does not change the conclusion, and none 
of the high-luminosity pulsars is affected. Details 
will be given elsewhere. 

The results shown in Fig. 1 demonstrate that a large 
fraction of high-luminosity pulsars correlate with 
sources of VHE 7-rays, emitting with a 7-ray lumi- 
nosity of order 1 % of the pulsar spin-down power 
The positive correlation does not necessarily im- 
ply that the pulsar or PWN itself is responsible 
for the 7-ray flux. It could also result from some 
other mechanism correlated with the pulsar or its 
creation, such as a supernova shock wave. The 
correlation found between 7-ray detectability and 
spin-down flux E/cP argues in favour of a pulsar- 
related origin of the 7-ray signal. On the other 
hand, for the PMPS pulsar sample, E/cP also cor- 



relates closely with the spin-down age T of the pul- 
ssii,E/cP r~^/^, and obviously with distance d, 
both parameters relevant for determining the 7-ray 
flux from shock-wave driven supernova remnants. 

The exact relation between pulsar parameters and 
7-ray luminosity is an interesting issue. Variations 
in exposure and hence in detection threshold over 
the survey range, as well as the uncertainty in pul- 
sar distance will smear out the turn-on curve of de- 
tectability versus E/cf shown in Fig. 1, but can- 
not fully account for the rather slow turn-on over 
a range of more than one order of magnitude in 
E/(P, combined with a detection probability be- 
low unity for even the highest-power pulsars. This 
indicates that E/cP cannot be the only parameter 
relevant for the 7-ray flux. The same conclusion 
is obtained from the observed variation of L^/E 
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of about an order of magnitude among the detected 
pulsars. However, this present pulsar sample is too 
small to investigate the dependence of on mul- 
tiple pulsar parameters, e.g. including pulsar age. 

A constant L^/E is also not necessarily expected. 
For a given age, the integral energy fed by the 
pulsar into the PWN increases with E. Apart 
from expansion losses, pulsar spin-down power is 
shared between particle energy and magnetic field 
energy. If equipartition between the two energy 
densities is assumed [13, 12], the magnetic field 
in the PWN will increase with E and hence the 
energy loss by synchrotron radiation will increase 
relative to and at the expense of inverse Comp- 
ton 7-ray production. Indeed, un-pulsed X-ray lu- 
minosity of pulsars is observed to increase faster 
then E, Lx oc £;i-4±o i [7] jjj ^ych scenarios, 
magnetic field values and therefore the balance be- 
tween X-ray and 7-ray emission will also depend 
on volume, i.e. on the expansion speed of the neb- 
ula and hence on the ambient medium. In addition, 
the current spin-down luminosity E may not be the 
only relevant scale; if the pulsar age is shorter than 
or comparable to the electron cooling time, relic 
electrons injected in early epochs with higher spin- 
down power will still contribute and may enhance 
significantly compared to the quasi-steady state 
achieved for old pulsars. 
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